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PREFACE 

In  a  previous  Memorandum,  RM— 288  7— PR,  "Minimum— Weight 
Design  for  Moving  Loads"  by  O.  Gross  and  W.  Prager,  a  problem 
involving  a  horizontal  I— beam  which  is  built  in  at  one  end  and  simply 
supported  at  the  other  was  proposed  and  solved.  The  rigorous 
mathematical  treatment  is  contained  in  RM— 2993— PR,  "A  Linear 
Program  of  Prager's:  Notes  on  Linear  Programming  and  Ex¬ 
tensions  —  Part  60."  The  present  Memorandum  considers  a  second 
beam— design  problem,  namely,  the  problem  of  a  beam  built  in  at 
both  ends  and  subject  to  the  same  loading  condition,  i.e.  ,  a  con¬ 
centrated  vertical  unit  load  moving  slowly  along  the  span. 

Beam  designs  of  the  type  (variable  cross  section)  considered 
in  this  Memorandum  are  perhaps  difficult  to  manufacture.  They 
nevertheless  represent,  in  a  sense,  the  ultimate  in  weight  saving 
and  are  therefore  at  least  of  theoretical  importance,  particularly 
in  the  view  of  possible  future  application  to  interplanetary  structures 
in  which  excess  weight  could  spiral  transportation  costs. 


V 


SUMMARY 


This  Memorandum  treats  the  problem  of  the  design  of  a  hori¬ 
zontal  I— beam  of  variable  flange  thickness  and  constant  web  height. 
The  beam  is  built  in  at  both  ends  and  is  otherwise  unsupported,  but 
is  subjected  to  a  concentrated  vertical  unit  load  that  moves  slowly 
from  one  end  to  the  other.  A  linear  program  in  function  space  is 
obtained,  using  the  concepts  of  limit  analysis  (based  on  the  signum— 
type  stress— strain  diagram),  and  the  solution  of  the  program  is 
given,  yielding  the  minimum  weight  design  as  well  as  the  critical 
reactions  at  the  ends  of  the  beam. 
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MINIMAL  WEIGHT  DESIGN  FOR  A  BUILT-IN  BEAM 


.  FORMULATION  OF  THE  PROBLEM 

Some  time  ago  W.  Prager  proposed  the  following  problem. 

An  I— beam  is  to  be  designed  to  fulfill  these  specifications:  The 
beam  is  to  be  built  in  at  both  ends  to  a  solid  structure.  The  web  is 
to  have  constant  height  and  thickness.  A  vertical  concentrated  load 
of  fixed  intensity  is  assumed  to  move  slowly  across  the  span.  As¬ 
suming  that  the  bending  moment  at  each  cross  section  is  taken  up 
by  the  flanges,  which  are  allowed  to  vary  in  cross-sectional  area, 
what  distribution  of  flange  cross  section  will  yield  a  minimal  weight 
design,  subject  to  the  constraint  that  the  beam  will  support  the  load 
in  plastic  flow? 

There  is  no  loss  of  generality  in  assuming  the  beam  to  be  of 
unit  length.  In  order  to  phrase  the  problem  in  mathematical  terms, 
it  is  helpful  to  have  recourse  to  the  following  Fig.  1,  depicting  the 
idealized  beam  together  with  its  variable  load  and  (statically  inde¬ 
terminate)  reactions.  We  assume,  also,  without  loss  in  generality, 
that  the  moving  load  is  of  unit  intensity. 


1 


1 


Fig.  I 
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Since  the  load  is  moving  sufficiently  slowly  to  obviate  dynami¬ 
cal  considerations,  the  beam  will  maintain  static  equilibrium. 
Imposing  the  force  and  moment  conditions  in  terms  of  the  load  and 
the  indeterminate  reactions  g^,  at  the  ends  (which,  for  a  fixed 
design,  are  functions  of  the  position  of  the  load)  then  gives 


(1. 

1) 

go(y>  +  gj(y)  =  1, 

(1. 

2) 

y  -  M^(y)  +  M^(y)  -  g^(y)  =  0 

(Moments  are  taken  counterclockwise  at  x  =  0.) 

Next,  we  determine  the  criticality  conditions.  If  the  critical 
design  bending  moment  of  the  beam  at  a  variable  position  x  is 
denoted  by  f(x),  then  in  order  that  the  beam  will  sustain  the  load 
in  an  arbitrary  position,  we  require  that  the  critical  bending 
moment  shall  nowhere  be  exceeded  in  absolute  value.  Two  cases 
present  themselves. 

Case  1.  0  <  y  <  X  <  1  (as  depicted  in  Fig.  1).  In  this  case, 

the  bending  moment  at  x  is  given  (looking  at  the  section  of  the 
beam  to  the  left  of  x  and  taking  moments  counterclockwise)  by 

M(x,  y)  =  X  -  y  —  xg^(y)  +  M^(y). 

Thus,  for  this  case,  the  criticality  condition  becomes 

(1.3)  f(x)  >  I  X  —  y  —  xg^(y)  +  M^(y)  |ifO<y<x<l, 


where  the  bars  denote  absolute  value. 
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Case  2.  0<x<y<l.  In  this  case,  upon  making  the  ap¬ 

propriate  modification  in  Fig.  1,  for  the  bending  moment  at  x 
(looking  at  the  section  of  the  beam  to  the  right  of  x  and  taking 
moments  clockwise)  we  obtain 

M(x,  y)  =  y-  x  —  (1-x)  gj(y)  +  M^(y). 

From  the  equilibrium  conditions  (1.1)  and  (1.2)  we  can  express  M 
in  terms  of  the  reactions  at  0;  thus,  a  little  algebra  gives 

M(x,  y)  =-xg^(y)  +  IVI^(y) 

(which  we  could  have  derived  anyway,  of  course,  on  sighting  to  the 
left  of  x).  Thus,  for  this  case,  the  criticality  condition  becomes 

(1.4)  f(x)  ^  i  -  xg^(y)  +  M^(y)  |  if0<x<y<l. 

Next,  observe  that  (1.  3)  and  (1.4)  can  be  combined  into  an 
equivalent  two— parameter  family  of  constraint  conditions,  namely 

(1.5)  f(x)  >  1  max(0,  x  -  y)  —  xg(y)  +  M(y)  1  0  5^5^- 

(We  have  dropped  the  subscript  0  for  convenience.) 

Finally,  since  the  web  of  the  idealized  beam  is  of  constant 
height,  the  design  bending  moment  at  x  will  be  proportional  to  the 
cross-sectional  area  of  the  flanges  there.  Hence,  the  variable 
portion  of  the  weight  of  the  beam  will  be  proportional  to 
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/  f(x)dx.  Since,  moreover,  it  follows  from  a  fundamental 

principle  of  limit  analysis  that  if  there  exist  reactions  M,  g  such 
that  the  bending  moment  is  nowhere  exceeded,  then  the  beam  will 
sustain  the  load  in  plastic  flow.  Thus,  the  mathematical  formulation 
of  the  design  problem  can  be  phrased  as  follows. 

Find  real  measurable  functions  f,  g,  M  on  [O,  l]  such  that 

.1 

/  f(x)  dx  =  minimum  , 

"  o 

subject  to 

(1-6)  f(x)  >  i  max(0,  x  —  y)  —  xg(y)  +  M(y)  | ,  0  5  y  £  ^  • 

The  solution  to  this  (continuous  linear  programming)  problem  will 
be  described  in  the  following  section.  Feasibility  and  minimality 
will  be  verified  in  Secs.  3  and  4,  respectively. 

2.  DESCRIPTION  OF  THE  SOLUTION 

We  now  assert  that  the  following  is  a  solution  to  the  design 
problem  as  formulated  at  the  end  of  Sec.  1: 


(2.  1)  f(x)  =  I  +  ^  log  V2(l  -  2.x)  if  0  <  X  <  i  , 


* 

Assuming  the  material  to  be  of  uniform  strength  and  density. 
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and  otherwise 


f(x)  =  f(l  —  x)  ; 


(2.2) 

and 

(2.3) 


g(y)  =  I  -  sinh  ^  if  f  1  y  5  I  ^ 


M(y)  =  I 


v 


if  T  <  y  <  V 

4  —  —  4 


The  value  of  the  program  is  given  by 
r  1  1 

(2.4)  j  f(x)dx  =  -  log  2  0.08664. 

o 

Graphs  of  these  functions  are  appended  at  the  end  of  this  paper. 

A  few  remarks  are  in  order  regarding  the  proposed  solution  and  its 
graphs . 

Remark  1 . 

We  observe  from  the  formulas  and  graphs  (particularly  Fig.  2) 
that  the  solution  is  symmetric.  This  is  not  surprising, since  it  can 
be  verified  a  priori  from  the  symmetry  of  the  model  and  linearity 
that  if  there  exists  a  solution  at  all,  there  exists  a  symmetric  one. 
Thus,  in  seeking  a  solution,  we  could  have  imposed  with  impunity 
the  condition 

Mo(y)  =  M^(l  -y) 

(upon  reverting  to  an  original  notation  of  Fig.  1).  Using  the  external 
equilibrium  conditions  (1.1)  and  (1.  2),  we  obtain 


-6- 


gjy)  =  1  -  y  +  M^(y)  -M^d  -y). 

Dropping  the  subscript  0,  as  before,  and  substituting  in  the 
constraint  condition  (1.6),  we  have 

f(x)  >  I  max(0,  X  —  y)  —  X  +  xy  +  xM(  1  —  y)  +  ( 1  —  x)  M(y)  |  , 

which  involves  only  one  unknown  reaction  function.  Thus,  if  x  >  y, 
we  must  have 

f(x)  >  I  —{1  —  x)y  +  xM(l  —  y)  +  (1  —  x)  M(y)  |  , 
and  if  x  <  y, 

f(x)  >  I  -(1  -y)x  +  xM(l  -y)  +  (1  -  x)  M(y)  |  . 

These  last  two  inequalities  are,  in  effect,  in  the  form  of  the  con¬ 
straint  conditions  in  the  mathematical  formulation  as  obtained  by 
Prager . 

Remark  2. 

Notice  that  we  have  not  defined  the  g  and  M  functions  outside 
the  interval  [1/4,  3/4]  (Figs.  3  and  4).  The  reason  for  this  is 
that  the  beam  does  not  attain  critical  condition  until  the  load  reaches 
the  middle  section  of  the  beam.  In  other  words,  the  reactions  are 
relatively  indeterminate  up  to  this  point;  i.e.  ,  infinitely  many 
feasible  extensions  of  these  functions  are  possible  outside  the 
critical  range.  We  shall  not  characterize  all  of  them  as  was  done 
in  [2]  ,  since  we  need  only  exhibit  one  such  extension,  as  is  re¬ 
quired  for  the  verification  of  feasibility  in  the  subsequent  section. 

As  it  turns  out,  simply  extending  the  functions  analytically  is  not 
feasible. 
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Remark  3. 


The  dotted  lines  in  Fig.  2  constitute  the  graph  of  the  bending 
moment  in  the  beam  (as  reflected  across  the  x— axis)  corresponding 
to  a  typical  position  of  the  load.  The  beam  is  critical  at  three  points 
(indicated  by  the  dots):  two  points  of  tangency  in  the  outer  sections, 
and  one  at  the  point  of  application  of  the  load.  The  graph  of  these 


critical  points  is  illustrated  in  Fig.  5.  The  two  curved  arcs  shown 


are  from  the  two  branches  of  the  hyperbola  y  =  x 


1 

4(2x-l)  ■ 


3.  VERIFICATION  OF  FEASIBILITY 

In  this  section  we  shall  give  a  partial  verification  that  with  an 
appropriate  extension  of  the  given  g  and  M  functions,  the  con¬ 
straints  (1.6)  are  met  by  the  solution  in  the  preceding  section.  To 
this  end,  we  extend  M  outside  [1/4,  3/4]  as  the  continuous  piece¬ 
wise  linear  extension  such  that  M(0)  =  M(l)  =  0.  This  extension  is 
shown  by  the  dotted  lines  in  Fig.  4.  Imposing  the  symmetry  con¬ 
dition 


g(y)  =  1  —  y  +  M(y)  -  M(  1  -  y) 

gives  rise  to  the  linear  extensions  of  g  shown  by  the  dotted  lines 
in  Fig.  3. 

We  are  now  required  to  verify  that  the  inequality 

f(x)  >  1  max(0,  X  —  y)  -  xg(y)  +  M(y)  1 ,  0  5  y  5 

is  satisfied  by  our  proposed  solution.  Our  first  observation  is  that 
because  of  the  symmetry  of  our  solution,  it  is  sufficient  to  verify 
this  inequality  for  ye  ( 0,  1/2]  .  Thus,  to  check  it  outside  the 
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interval  ^1/4,  3/4^  ,  we  need  only  consider  the  case  0  <  y  <  1/4. 
For  this  case  however,  and  for  any  given  value  of  x,  the  function 
of  y  under  the  absolute  value  sign  is,  by  our  extension,  piecewise 
linear  in  at  most  two  pieces,  and  hence  will  attain  its  maximum 
absolute  value  at  one  of  the  points  y=0,  y=x  or  y  =  1/4.  Since 
g(0)  =  1  and  M(0)  =  0,  the  case  y  =  0  reduces  to  f(x)  >  0.  We 
leave  the  straightforward  verification  of  this  inequality  as  an 
exercise.  Next,  if  y  =  x,  then  in  this  case  we  must  have 
0  <  X  <  1/4,  and  hence  we  are  required  only  to  verify  that 

f(x)  >  1  — xg(x)  +  M(x)  1  if  0  <  X  <  1/4. 

But  one  checks  for  this  case  that  the  function  on  the  right  is  simply 
a  quadratic  function  which  vanishes  at  x  =  0  and  x  =  1  suid  hence 
will  attain  its  maximum  on  the  specified  range  at  x  =  1/4/  On  the 
other  hand,  since  f  is  decreasing  on  this  interval,  in  order  to 
check  the  inequality  it  suffices  to  check  it  at  x  =  1/4;  i.e.  ,  it 
suffices  to  verify  that 

>|-I6  (j)+M(i)|  . 

However,  using  formulas  (2.1),  (2.2),  and  (2.3),  which  now  apply, 
readily  gives 

<l)  '  |-ie(j)  +m(7)|- 

This  disposes  of  the  case  y  =  x  <  • 


-9- 


Finally,  we  note  that  the  case  y  =  1/4  is  included  by  continuity 
in  the  range  [1/4,  3/4]  .  Thus,  to  complete  the  verification  it  is 
sufficient  to  establish  that 

f(x)  >  I  max(0,  X  —  y)  —  xg(y)  +  M(y)  |  if  0  <  x  <  1 

and  1/4  <  y  <  3/4.  In  this  range,  however,  formulas  (2.2)  and 
(2.  3)  apply.  Now  observe,  that  for  a  given  value  of  y  in  this 
range,  the  function  of  x  under  the  absolute  value  sign  is  piecewise 
linear  (in  at  most  two  pieces).  We  assert  that  f  is  simply  the 
absolute  upper  envelope  of  this  family  of  functions,  i.  e.  , 


f(x)  =  max  |max(0,  x  —  y)  —  xg(y)  +  M(y)  |  0  <  x  <  1. 


Since  the  evaluation  of  the  right  member  of  the  above  was  precisely 
the  method  by  which  the  function  f  was  obtained,  we  shall  leave 
this  elementary  albeit  arduous  calculus  exercise  to  the  reader, 
lest  this  section  become  as  long  as  the  corresponding  one  in  [2]  . 


4.  VERIFICATION  OF  MINIMALITY 


In  this  final  section  we  shall  show  th  ,1  the  proposed  feasible 

solution  given  by  (2.  1) — (2.  3)  satisfies  the  minimality  requirement 

and  hence  solves  the  problem  proposed  in  Sec.  1. 

.1 


Let 


1_ 

4 


dx 


that 


A 

4 


2x  -  1 


It  is  then  a  trivial  matter  to  check 


J*  log  2. 
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Next,  by  means  of  well  known  identities  on  the  inverse  hyperbolic 
sine,  it  becomes  a  straightforward  exercise  to  verify  that  the 

* 

functions  f,  g,  M  given  in  Sec.  2  satisfy  the  following  relations: 


(4.  1) 


t(x)  +m(x-.5^)  it  0  <  X  <  i 


=  xg(x)  —  M(x) 


li  -  <  X  <  — 
4  —  —  4 


Upon  making  the  changes  of  variable  yielding  x*  =  x  —  l/4(2x  —  1) 
in  the  1st  and  3rd  integrals  in  the  right  member  of  the  above, we 
obtain,  respectively  (using  the  appropriate  branches  of  the  hyper¬ 
bola) 

* 

Notice  that  the  values  of  the  arguments  of  the  g  and  M 
functions  used  lie  in  the  appropriate  range  (1/4,  3/4). 
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1 


Using  the  fact  that  the  two  outer  integrals  are  now  expressed  as 
integrals  over  the  same  interval  ^1/4,  3/4^  ,we  can  add  their 
integrands,  collecting  coefficients  of  g(x')  and  M(x‘),  to  obtain 
as  their  joint  contribution: 


3 


(— x’g(x')  +  M(x'))  dx’  . 


4 
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Upon  substituting  this  joint  contribution  in  (4.  2)  we  obtain 


14  4 

(4.3)  r  f(x)dx  =  J*  +  r  (xg(x)  —  M(x))dx  +  j  (— x'g(x')  +  M(x'))dx' 
^o  '^l  ^1 


since  the  integrals  involving  the  reactions  cancel. 

Next,  let  (f  ,  g,  M)  denote  any  feasible  triple  of  functions 
satisfying  the  constraints  (1. 6),  i.e.  , 


(4.  2)  f  (x)  >  lmax(0,  x  —  y)  —  xg(y)  +  M(y)  1  ,  0  <  <1 

y 


In  particular,  we  must  have 


(4.3)  f(x)  >  max(0,  x  —  y)  —  xg(y)  +  M(y),  0  <  £  1- 

y 

Now,  one  verifies  that  on  the  interval  (0,  1/4) 


(4.4)  0  1^1  ^ - 2x  -  1  - 


Hence  (4.  3)  must  hold  for  y  =  x  -  i.e.,  using  (4.  4),  we  must 

have 


(4.5)  ?(x)  >  -  xg(x  -  2x1  1  )  +  M  (x  -  ^  )  if  0  <  X  < 


rH 


-1 3- 


Next,  (4.  2)  implies 

f(x)  >-  max(0,  X  -  y)  +  xg(y)  —  M(y),  0  <  ^  <  1 

which  must  hold  for  y  =  x.  In  particular,  on  the  interval 
(1/4,  3/4),  we  must  have 

(4 . 6)  f(x)  >  xg(x)  —  M(x)  if  ^  _5  ^  ^  • 


Now  observe  that  (4.5),  (4.6),  and  (4.8)  are  precisely  the  relations 
(4.  1)  with  "  replaced  by  Consequently,  of  one  integrates 

the  inequalities  over  their  ranges,  adds  them  together,  and  makes 
the  same  changes  of  variable  as  before,  the  integrals  involving  the 
unspecified  reactions  g,  M  will  cancel  to  obtain 

f(x)dx  ^  J*  -  [  f(x)dx. 

^o 

Thus  the  design  is  indeed  minimal  and  the  verification  is  complete. 
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Perhaps  a  word  or  two  is  in  order  regarding  the  almost  com¬ 
pletely  synthetic  treatment  of  beam— design  problems,  since  the 
reader  is  no  doubt  wondering  how  the  author  came  to  obtain  the 
obviously  complicated  solution  in  the  first  place.  The  author  need 
not  apologize  for  this  insofar  as  he  has  indicated  a  complete  veri¬ 
fication.  Suffice  it  to  say,  however,  that  the  "method”  of  derivation 
involved  quite  a  bit  of  judicious  guessing,  solving  of  differential 
equations,  matching  end  conditions,  etc.  Lest  this  Memorandum 
become  unduly  long,  we  shall  not  go  into  details  here,  but  rather 
refer  the  reader  to  [l]  ,  where  motivation  for  solving  a  similar 
beam— design  problem  is  given. 

We  shall  close  this  Memorandum  with  some  remarks  about  the 


loads  that  our  beam  design  will  sustain.  Let  W  denote  an  arbitrary 
downward  static— load  distribution,  of  total  weight  not  exceeding 

f  1 

unity,  i.e.,  dW(y)  >  0  and  j  dW  (y)  <  1.  We  assert  that  our 

beam  will  support  this  load.  This  is  perhaps  intuitively  clear. 
Nevertheless,  we  shall  present  a  rigorous  verification.  The 
appropriate  modification  of  Fig.  1  yields  the  following  as  the 
criticality  condition  corresponding  to  (1.6); 


(4.9)  f(x)>|  [  (x-y)dW(y)  -xgj^(W)  +  M^(W)|,  0<x<l, 

-  o 


where  g^  and  M  (W)  are  the  support  and  moment  reactions, 

Lj(W)  l# 

induced  by  the  load,  at  the  left  end  of  the  beam.  Thus,  the  beam 
will  sustain  W  if  there  exist  nonnegative  functionals  g^  and  Mj^ 
satisfying  (4.9).  Now,  we  have  already  verified  that 


f(x)  >  max(0,  x  —  y)  —  xg(y)  +  M(y),  0  £  ^ 
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Integrating  this  last  (over  [O,  l]  )  with  respect  to  dW(y)  >  0  and 

noting  f(x)  _>  0  and  /  dW(y)  <  1,  we  obtain 

'-'o 


f 


X 


f(x)  >  (x  — y)dW(y)  — X  /  g(y)dW(y)+  /  M(y)dW(y). 


J. 


"o 


Similarly,  using  the  established  inequality 


f(x)  >  —  max(0,  X  —  y)  +  xg(y)  —  M(y), 

we  obtain 

X  1  ,1 

f(x)  >  —  /  (x  —  y)  dVv  (y)  +  X  /  g{y)dW(y)—  j  M(y)dW(y), 

Thus,  it  suffices  to  set 

g  (W)  =  '  g(y)d\V(y) 

^  o 

and 

M  (W)  =  [  M(y)dW(y). 

^  ^o 

Combining  the  two  inequalities  then  gives  (4.9).  Hence  the  beam 
will  sustain  the  load.  Moreover,  since  criticality  is  not  reached, 
for  examiple,  by  a  unit  concentrated  load  until  it  reaches  its  mid 
section,  the  beam  will  support  other  loads  as  well.  We  shall 
refrain,  however,  from  giving  a  complete  characterization.  Thus, 
in  effect,  we  have  solved  a  more  general  problem  than  the  one 
proposed!  We  can,  moreover,  phrase  our  result  in  another  way. 
Let  r  be  a  family  of  downward  load  distributions  over  [O,  l]  , 
i.  e.  ,  W  €  dW(y)  >  0  if  Then  there  exists  a 
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built— in  beam  design  capable  of  supporting  any  W  e  such  that 
its  design  bending  moment  (}>  satisfies 

<j)(x)  dx  <  m  sup.  /  dW(y), 

W  € 

where  m  =  -^  log  2  (provided  the  right  member  is  finite),  the 

O 

constant  m  being  best  possible. 

A  curious  research  problem  suggests  itself  at  this  point: 

Given  a  built-in  beam  design  B  of  specified  variable  cross  section, 
find  a  beam  design  of  minimal  weight  which  will  sustain  every  down¬ 
ward  static— load  distribution  that  B  does.  It  is  not  immediately 
clear,  for  example,  that  B  is  itself  minimal.  We  note  in  passing, 
however,  that  for  the  case  of  a  beam  that  is  simply  supported  at 
each  end  and  of  concave  design  moment,  the  assertion  is  true. 

As  an  amusing  sidelight  on  the  foregoing  remarks,  consider 
the  situation  in  which  two  men  of  weight  one— iialf  unit  each 
(Martians,  no  doubt)  are  walking  across  our  beam  in  opposite 
directions.  Since  the  beam  will  not  become  critical  unless  they 
happen  to  pass  each  other  in  its  midriff,  the  little  fellows  should 
be  careful  to  be  out  of  step  in  case  of  such  an  event. 

* 

More  generally,  if  its  (continuous)  graph  is  "star  shaped" 
with  respect  to  the  points  (0,  0)  and  (1,  0),  then  this  last  con¬ 
dition  is  both  necessary  and  sufficient. 
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Design  moment  of  beam 


Vertical  reaction  at  left  end  of  beam 


Criticol  points  in  beam  (plastic  hinges  form  at  these  points) 


X 


Position  of  load 


Fig.  5 — Critical  points  in  beam  versus  load  position 
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